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LONG-TERM  GOAL 

The  lack  of  good  ocean  snrface  data  on  the  one-mile  scale  is  a  major  stnmbling  block  for  proper 
evalnation  of  proposed  Mobile  Offshore  Base  (MOB)  conceptnal  designs.  A  general  model  and 
associated  software  for  stationary,  linear,  statistical,  directional  wave  systems  have  been 
available  for  some  time.  If  that  model  can  appropriately  be  nsed  for  design  of  marine 
mega-strnctnres,  then  there  is  no  problem  in  proceeding  with  conceptnal  MOB  designs  and  their 
evalnations.  However,  it  is  known  that  storm  waves  are  nonlinear  to  some  extent.  The  major 
goal  is  the  determination  of  the  accnracy  of  the  linear  approximation  and  the  type  of  nonlinear 
featnre  present  in  severe  storm  waves. 

OBJECTIVES 

Many  frnstrations  and  nncertainties  wonld  disappear  if  adeqnate  data  on  the  proper  spatial 
scale  were  available  to  serve  as  a  basis  for  selection  of  design  criteria.  This  project  is  designed  to 
examine  and  analyze  NASA  scanning  radar  altimeter  data  collected  by  Ed  Walsh. 

APPROACH 

Walsh  has  recorded  SRA  data  in  flights  over  storm  wave  systems  near  Anstralia  and  is  in  the 
process  of  collecting  data  for  Atlantic  storms  (1998  and  1999  hnrricane  seasons).  In  the  present 
project,  these  data  are  being  analyzed  and  snmmarized  graphically  in  a  “design  engineer’s  atlas 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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of  storm  wave  measurements” .  The  atlas  includes  various  statistical  studies  that  address  the 
linear/nonlinear  balance  in  severe  storm  waves. 

WORK  COMPLETED 

The  work  this  year  consisted  of  developing  methodology  and  algorithms,  and  related  software. 
Pre-existing  storm  measurements  and  data  from  several  Atlantic/Gulf  of  Mexico  hurricanes 
were  analyzed. 

RESULTS 

The  work  leading  up  to  the  VLFS’99  meeting  in  September  1999  encountered  questions  which 
needed  resolution  before  the  production  of  the  final  statistical  Atlas  of  Storm  Topography.  In 
particular,  methods  were  needed  which  did  not  interpolate  to  a  regular  grid  and  which  allowed 
computations  to  be  based  on  the  actual  measurements  no  matter  how  they  were  irregularly 
spaced  over  [x,  y).  This  was  important  in  order  to  guarantee  that  no  local  nonlinear  features 
actually  present  in  the  sea  surface  were  seriously  reduced  or  even  eliminated  in  the  data 
processing.  Some  of  the  procedures  previously  used  varied  somewhat  with  the  orientation  of  the 
regular  grid  which  was  tied  to  the  flight  direction.  Consequently,  this  is  also  an  artifact 
introduced  by  the  analysis  method,  so  it  was  desirable  to  find  procedures  that  were  insensitive 
to  flight  direction  and  which  would  work  even  if  the  airplane  were  Hying  “S”  patterns  or  sharp 
turns  within  the  data  section  being  analyzed. 

Two  methods  were  selected  for  serious  review.  These  were  distorted  quadrilaterals  and 
Delaunay/  Voronoi  tesselations  (Ahuja  and  Schachter,  1983).  The  Delaunay /Voronoi  approach 
was  judged  to  have  the  most  generality  for  data  with  irregular  spacing,  and  conveniently,  could 
be  computed  with  functions  DELAUNAY  and  VORONOI  currently  existing  in  MATLAB  5.3. 

A  tesselation  is  any  family  of  disjoint  2-D  patterns  covering  exhaustively  a  selected  planar  set. 
Let  {x^,  y{,  i  =  1,2, ,  AI}  be  a  set  of  locations  on  the  plane.  The  ^th  Voronoi  set  is  the  set  of  all 
points  in  the  plane  which  are  closer  to  the  (a;^,  y^j-location  than  to  any  of  the  other  specified 
locations.  This  turns  out  to  be  the  polygon  formed  by  the  perpendicular  bisectors  of  the  line 
segments  from  (a;^,  to  the  other  neighboring  locations.  An  example  of  Voronoi  polygons  is 
given  in  Fig.  1.  The  red  plus  symbols  are  the  generation  locations,  and  each  polygon  surrounds 
those  points  which  are  closer  to  the  enclosed  location  than  to  any  other  of  the  locations.  The 
area  of  the  Voronoi  polygon  is  the  natural  value  to  assign  to  dx  dy  in  the  numerical  integration 

^  =  f  i^v  Vi)  dx  dy 

over  irregularly-located  measurements.  It  allows  a  reasonable  definition  of  a  two  dimensional 
finite  Fourier  transform  over  irregularly-spaced  data. 

The  Delaunay  triangulation  consists  of  the  triangles  formed  by  the  line  segments  perpendicular 
to  the  sides  of  the  Voronoi  polygons  connecting  the  two  locations  which  generated  that  side. 

The  Delaunay  triangles  corresponding  to  the  Voronoi  polygons  in  Fig.  1  are  shown  in  Fig.  2. 
[Please  note:  The  Delaunay  line  segment  may  not  physically  intersect  the  polygon  side  which  it 
is  perpendicular  to.) 


Each  of  the  locations  not  on  the  convex  hnll  of  the  set  of  locations  is  completely  snrronnded  by 
Delannay  triangles  and  a  set  of  radiating  line  segments  which  form  the  sides  of  the  snrronnding 
triangles  that  are  connected  to  that  location.  The  snrronnding  set  of  verticies  connected  to  a 
specified  location  by  a  single  Delannay  line  segment  will  be  called  the  D-neighborhood  of  order 
one,  or  primary  D-neighborhood,  for  the  selected  locations.  Similarly,  the  set  of  verticies  not  in 
the  primary  D-neighborhood  and  connected  by  two  Delannay  line  segments  to  the  specified 
location  will  be  called  its  secondary  D-neighborhood,  or  D-neighborhood  of  order  two.  This 
concept  can  be  extended  to  neighborhoods  of  order  3,  4,  etc.  in  an  expanding  set  of  neighboring 
locations  snrronnding  the  reference  locations. 

There  is  a  growing  body  of  mathematical  techniqnes  bnilt  aronnd  nearest-neighbor  concepts. 
Many  of  these  are  based  on  Enclidean  distance  or  some  other  measnre  of  distance,  and  the 
neighborhood  of  locations  often  end  np  all  on  one  side  of  the  reference  location.  The 
D-neighborhoods  avoid  this  since  they  always  snrronnd  the  reference  location  if  it  is  not  on  the 
convex  hnll. 

In  the  stndy  of  the  irregnlarly-spaced  SRA  data,  the  D-neighborhoods  form  a  natnral  setting  for 
the  stndy  of  local  nonlinearities.  The  elevations  of  each  location  are  actnal  SRA  measnrements, 
so  there  is  no  distortion  introdnced  by  interpolation. 

The  Delannay  triangles  can  be  nsed  with  splines  to  define  the  snrface  more  precisely.  In  the 
work  over  the  last  several  months  cnbic  qnadratic,  and  planar  splines  were  examined  as  possible 
choices  to  nse  with  the  finite  elements  made  np  of  the  Delannay  triangles.  This  leads  in  a 
natnral  way  to  contonring  rontines  for  elevation  based  on  the  triangnlation,  and  also  contonr 
maps  of  snch  sea  snrface  properties  as  slope. 

Since  qnadratic  and  cnbic  splines  may  introdnce  artifacts  of  analysis  difficnlt  to  evalnate,  it  was 
decided  to  stick  with  the  most  simple  case  of  planar  finite  elements,  at  least  for  the  present. 
Thns  on  each  Delannay  triangle,  with  vertices  i  =  1,2,3},  the  eqnation  for  the  elevation 

is 

z  —  A  Bx  -\-  Cy 

with 

Zi  —  A  Bx\  A'  Cy\ 

^2  —  A  T  Bx2  -\~  Cy<2 
Z3  —  A  A-  Bx3  A-  Cy^ 

or 

Zi  1  \  I  Xi  y-y^\  A 

Z2  ^  1  X2  y2  B 
Z3  \  [1  X3  y^  \  [  C 

The  constants  A,  B,  and  C  for  that  finite  element  may  be  solved  for  with 

1  1-1  r 

1  xi  y^  zi 

1  2:2  y2  Z2 

1  %  J  L  ^3 


providing  the  triangle  is  not  degenerate. 


Let  (y4i,  5i,  Cl)  and  {A2,  B2,  C2)  be  the  constants  for  two  adjacent  Delannay  triangles.  Then  the 
angle  between  them  is  given  by 


B1B2  +  Cl  C2 

cos  w  =  — I 

+  c?)  (i^  +  c5) 

A  contonr  map  of  cos  9  shonld  delineate  crests  and  tronghs  fairly  natnrally,  since  this  is  where 
the  change  in  slope  is  largest. 

Another  choice  for  locating  the  crest  and  trongh  lines  is  based  on  the  Delannay  line  segments. 
Each  segment  is  common  to  two  triangles.  The  slope  of  z  perpendicnlar  to  the  segment  in  each 
of  the  two  triangles  can  be  nsed  to  compnte  a  directional  change  in  slope  perpendicnlar  to  the 
line  segment.  A  line  segment  can  be  colored  in  a  MATLAB  fignre  according  to  the  change  of 
slope  across  that  segment.  A  crest  line  wonld  then  be  made  np  of  the  Delannay  segments  with 
change  from  positive  slope  to  negative  slope  across  it.  This  appears  very  promising  in  the 
compntations  so  far. 

The  two  MATLAB  fnnctions,  DELAUNAY  and  VORONOI,  do  not  have  snpporting  fnnctions  of 
the  type  needed  in  the  SRA  wave  analysis.  Conseqnently  considerable  time  dnring  the  last  few 
months  have  been  devoted  to  the  design  and  programming  of  a  snite  of  fnnctions  for  the 
operations  on  the  Delannay  triangnlation. 

MATLAB  code  allows  nested  do-loops,  bnt  their  nse  with  large  array  compntations  eats  np 
mnch  compnter  time.  Conseqnently,  efficient  processing  of  large  data  sets  reqnires  programs 
containing  no  do-loops,  if  at  all  possible.  Sometimes  it  can’t  be  avoided.  However,  nsnally  with 
snfficient  insight  and  ingennity,  do-loops  can  be  avoided  by  working  with  matrix  and  array 
calcnlations.  The  work  on  this  snite  of  compntational  fnnctions  is  cnrrently  nnderway. 

IMPACT  /APPLICATION 

To  nnderstand  the  very  nonlinear  rogne  waves  we  mnst  first  record  the  conditions  and  patterns 
that  prodnce  them.  How  to  obtain  snch  data  was  a  topic  of  considerable  discnssion  at  the 
September,  1999  VLFS  meeting.  These  and  snbseqnent  discnssions  indicate  that  we  shonld 
develop  a  mechanism  for  monitoring  and  recording  form  and  conditions  of  the  air-water 
interface  in  three  dimensions  and  over  a  relatively  long  period  of  time.  To  this  end,  we  have 
proposed  a  snpplemental  project  (Marrs  et  al,  2000)  that  will  attempt  to  obtain  rapid 
seqnences  of  stereographic  images  from  which  sea-snrface  topographic  data  can  be  extracted  via 
digital  photogrammetry.  The  resnlting  data  sets  will  provide  3-D  representations  of  waves 
throngh  a  specified  time  period.  We  propose  to  record  and  analyze  conditions  prior  to,  dnring, 
and  following  the  departnre  of  wave  patterns  from  linear  behavior. 
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